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Classical and quantum dynamics of a harmonic oscillator in a monochro- 
matic wave is studied in the exact resonance and near resonance cases. This 
model describes, in particular, a dynamics of a cold ion trapped in a linear 
ion trap and interacting with two lasers fields with close frequencies. Ana- 
lytically and numerically a stability of the "classical ground state" (CGS) - 
the vicinity of the point (x = 0,p = 0) - is analyzed. In the quantum case, 
the method for studying a stability of the quantum ground state (QGS) is 
suggested, based on the quasienergy representation. The dynamics depends 
on four parameters: the detuning from the resonance, 5 = ^ — Vi/u), where O 
and Lo are, respectively, the wave and the oscillator's frequencies; the positive 
integer (resonance) number, the dimensionless Planck constant, h, and the 
dimensionless wave amplitude, e. For 5 = 0, the CGS and the QGS are un- 
stable for resonance numbers £ = 1, 2. For small e, the QGS becomes more 
stable with increasing 6 and decreasing h. When e increases, the influence of 
chaos on the stability of the QGS is analyzed for different parameters of the 
model, i, 6 and h. 
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I. INTRODUCTION 



One of tlie major difficulties in developing quantum technologies, such as quantum 
computers,0'l are different kinds of specifically quantum dynamical instabilities that can 
occur due to interactions between different degrees of freedom and resonant interaction with 
the external fields. Instabilities in quantum systems have different nature than instabilities 
in classical systems, in which such kind of phenomenon as dynamical chaos occurs as a result 
of exponential divergence of initially close trajectories. In quantum systems, the notion of 
a trajectory is not well defined. This is one of the main reasons why most of well-developed 
methods for stability analysis can not be directly applied to quantum systems. 

In this paper the dynamics of a harmonic oscillator in a monochromatic wave field is 
studied. This system describes, in particular, a quantum dynamics of a cold ion trapped in 
a linear ion trap and interacting with two laser fields with close frequencies.i Our attention 
is focused mainly on classical and quantum behavior in the region of parameters close to the 
quantum ground state (QGS) of the harmonic oscillator. This case corresponds to classical 
dynamics in the vicinity of the point (x = 0, p = 0) in the phase space. We shall call 
the corresponding region of parameters for the classical system a "classical ground state" 
(CGS) of a harmonic oscillator. The stability of the CGS and the QGS at stable and 
chaotic classical regimes is analyzed for different parameters of the model. In Sec. II, we 
consider classical dynamics in the vicinity of the CGS in the exact resonance case, 6 = 0. 
In this situation the system we study is degenerate,^ and an infinitely small perturbation 
generates in the classical phase space an infinite number of the "resonance cells" separated 
by the infinite stochastic web. A classical dynamics inside the resonance cells at small e 
is described by using the resonance perturbation theory.@ It is shown, that in the vicinity 
of the CGS (the "central cell"), and for small enough e, the dynamics in the case 5 = is 
unstable for resonance numbers £=1,2 and stable for i > 3. 

We show that in spite of the "resonant Hamiltonian" describes many features of the 
dynamics inside the resonance cells, it can not be used for describing the stability of the 
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system in the central cell, at i > 3. The classical dynamics in this region is determined, by 
the Mathieu equation. It is shown that at small enough value of e the area of the central 
cell increases with increasing the resonance number, i, and increasing the perturbation 
amplitude, e. 

The dynamics of the system near the CGS in the near resonance case, when 5 7^ 0, is 
considered in Sec. III. It is shown that in this case at small e the classical dynamics near 
the CGS is stable at any value of the resonance number i. The cases i = 1 and i = 2 are 
considered in detail. It is shown that the CGS in these two cases becomes unstable at much 
less values of e, than for large £ (£ = 4, 5...). 

A stability of the quantum system is considered in Sec. 4. In this case, an additional 
parameter, a dimensionless Planck constant, h, influences significantly the behavior of the 
system. Because the Hamiltonian is time-periodic, with the period 2tt/Q, we use the Floquet 
theory to study localization properties of quantum system in the region of the QGS. For 
e <^ 1, and small enough h, the quasienergy (QE) states can be conventionally divided 
into the following groups. The first group includes the QE states which belong to some 
particular "resonance cells" in the Hilbert space. These QE states are well-localized inside 
a given resonance cell. The second group includes the delocalized "separatrix" QE states 
which correspond to the classical stochastic web. Usually, these "separtrix" QE states are 
responsible for tunneling effects between different "resonance cells" , even for small e.i 

We show that the QGS is stable when the following conditions are satisfied: (a) an exis- 
tence of the QE state mainly localized in the QGS of the harmonic oscillator, (b) when e <^ 1 
and chaos is weak, (c) small enough h, when one can neglect the tunneling phenomenon. 
When h is larger than the size of the central cell, no QE state localized in the QGS of the 
harmonic oscillator was found. We show, that at small enough values of e, the stability 
of the QGS can be improved by choosing the non-resonant frequency of the wave, so that 

For small enough e, a stability of the QGS is mainly determined by the structure of the 
QE state mostly localized at the QGS of the harmonic oscillator (QGS QE state). In the 
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exact resonance case this particular QE state has a complicated structure. On the one hand, 
it is mostly localized on the QGS, but, on the other hand, it has the "separatrix" structurelll 
and provides a tunneling from the region of the QGS to other resonance cells. In order 
to improve stability of the QGS, the separatrix structure should be destroyed by choosing 

In the exact and near resonance cases we study numerically the probability, Pq, for a 
particle to remain in the QGS, depending on e and for different values of h, C. and 5. When 
(5 = and e is small, the dynamics depends on the value of the quantum parameter, h, which 
controls degree of delocalization of the QGS QE state. When h is small, the probability of 
tunneling to other cells is small too, and Pq decreases with e increasing, because chaos makes 
the QGS QE state more delocalized. For large enough values of /i, the dependence of Pq 
on e becomes more complicated because in this case we have two superimposing effects: the 
influence of chaos on the dynamics, and the effect of tunneling to the classically unacceptable 
cells. These two features make the dependence Pq = -Po(e) non-monotonic. In the region 
of small e, Pq increases with e increasing. This behavior of Pq can be explained as a result 
of destroying the delocalized QE states with increasing e. Further increase of e leads to 
decreasing Pq. This behavior of Pq is connected with delocalization of QE states due to 
increasing of the chaotic component. At large enough values of e, the dependence of Pq on 
e becomes more complicated, and generally must be considered using many QE states. It 
is shown that the value of S influence significantly the stability of the QGS only at small 
e < 1. The stability of the QGS at £ = 1, 2 and 5 7^ is explored. It is shown that the 
QGS becomes unstable at much less values of the wave amplitude than in the case when 
i is large. The derived in this paper results are important for understanding a stability of 
quantum dynamics in the vicinity of the ground state. 
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II. CLASSICAL DYNAMICS NEAR THE CGS IN THE CASE OF EXACT 

RESONANCE 



The Hamiltonian of the harmonic oscillator in a monochromatic wave is, 

H=^ + + vo cos{kx - Qt) = Ho + V{x, t), (1) 

where M is the mass of the particle, p is the momentum, k is the wave vector, vq is the 
amplitude of the perturbation, and Hq is the Hamiltonian of the harmonic oscillator. In the 
first part of this section we discuss the case of the resonance, when Q = iu. It is known 
(see for example Ref.0) that under the resonance condition the infinitely small perturbation, 
Vq, is enough to generate in the classical phase space the infinite stochastic web. The web 
is inhomogeneous, and its width decays with decreasing the perturbation amplitude, Vq, 
and increasing \p\ and Inside the cells of the web a particle moves along stable closed 
trajectories. 

To analyze the dynamics of the harmonic oscillator in a monochromatic wave, described 
by the Hamiltonian (|I]), it is convenient to use the resonance perturbation theory! discussed 
below. Let us perform a transformation from the variables {p,x) to the canonically conju- 
gated variables {J^, (p), 

X = {2 / M Luy/^ sin (f = r{J^) sincp, (2) 

Px = {2J^Mujy^'^ cosif = Mujr{J^) cos (f, (3) 

where r{J^p) = {2J^/ MuiY^'^ is the amplitude of oscillations. It is more convenient to work 
with the dimensionless coordinate, X = kx, and the dimensionless momentum, P = kp/Mu, 
which are related to the variables {J,^, ip) by the formulas, 

X = p{J^) sinif, (4) 
P = p{J^) cos ip, (5) 
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where p(J^) = \JX'^ + = kr[J^p). In order to treat the time on the same ground as the 
phase ip let us introduce the new pair of canonically conjugated variables, (J/3, where 
(3 = Qt. The initial Hamiltonian (|I]) expressed through the new variables takes the form, 

H = J^uj + JpQ + t>o cos (p sin ip — P) . (6) 

It is independent of time, but describes the motion in the two-dimensional space. The 
nonlinear perturbation in Eq. can be expressed in a series, 

00 

Vq cos (p sin p- p) =vo ^ J„(p) cos {nup - j3) , (7) 

m=— 00 

where Jm{p) is the Bessel function. Under the resonance condition, Q = iu, all terms in the 
sum d^) quickly oscillate and can be averaged out, except for one term with m = i. In this 
approximation, the Hamiltonian (|^) is reduced to, 

H = J^uj + JpVt + VqJi,{p) cos {Ip - (3) . (8) 

It is convenient to introduce new, resonance, variables, (/, 9), (J, by using the generating 
function, 

F = i{ip-p) + jp, 

The new Hamiltonian, 

H = i{iuj + jiu + voMp) cos e, (9) 

where 6 = ip — (3, is independent of the variable f3. Hence, J = const. The resonance 
Hamiltonian, 

Heip, e)=H - Juj = VoJeip) cos 9, (10) 

(where we used the resonance condition iu = Q) is independent of time, unlike the initial 
Hamiltonian (|1]). 

The Poincare surfaces of section of the system described by the Hamiltonian ([^) in 
variables {X, P) are shown in Figs. 1 (a) - 1 (e), for the cases i = 1,2, 3, 4, 5. The phase 
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FIG. 1. The resonance cells in the phase space for e = 0.05, 6 = and: (a) £ = 1, {h) i = 2, (c) 
e = S,{d)e = 4, (e) £ = 5. 
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points are plotted at the moments tj = 27ij/Q, where j = 0, 1, 2, . . .. One can see that 
the phase space has an axial symmetry of the order i. The phase space is divided into 
the cells. A particle moves along closed trajectories inside the cells. (In Figs. 1 (a) - 1 
(e) only the boundaries of the cells are demonstrated.) At small values of vq, the motion 
inside the resonant cells, illustrated in Figs. 1 (a) - 1 (e), can be considered in the resonance 
approximation. The next order approximation is required only for analysis of the motion 
inside the exponentially small chaotic regions near the separatrices. It is shown below that 
the resonance approximation also fails to describe the motion in the region near the point 
(X = 0, P = 0). 

The symmetry of the phase space in Figs. 1 (a) - 1 (e) follows from the form of the 
resonance Hamiltonian ([10|), which is invariant under the transformation, 

ip-^ip + 27c/i. (11) 

In the phase space there are i cells connected by the transformation (jl^), and i cells con- 
nected by the same transformation but described by the Hamiltonian Hilp, 6 + it) in Eq. 
(PUp. The total number of cells at given p is 2i. Thus, in Fig. 1 (a) the upper cell cor- 
responds to the Hamiltonian Hi{p,6), while the lower cell corresponds to the Hamiltonian 
Hi{p, 9 + n) in Eq. (|lO|) ; the upper and lower cells in Fig. 1 (b) correspond to i^2(p, 9), and 
the right and the left cells correspond to i^2(p, 9 + n), and so on. 



It is easy to see, the resonance Hamiltonian (10) yields unstable solution near the CGS 



(the point {X = 0, P = 0)). To show this, we present the Hamiltonian (10) in the form 



P 

He = vo^^ cosiip = El, (12) 

where Ei = const (because the resonance Hamiltonian is independent of time). Also, we 
took into account that at p ^ 1 the Bessel function can be expressed in the form, 

Mp) ^ p'/2'i\, (13) 

and the fact that in the Poincare surfaces of section the position of the particle is taken 



at the moments ^Itk = 27r/c, /c = 0, 1, 2, . . .. It follows from Eq. ([12D that at the angles 
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ifk = ||(2A; — 1), = 1, 2, . . . , 2£, the radius p should sharply increase or decrease. It is seen 
from Figs. 1 (a) - 1 (e) that at angles (pk the particle moves in the radial direction. The 
grows of p is limited due to nonlinear effects. 

However, the resonance perturbation theory does not adequately describe the motion of 
a particle in the vicinity of the point (X = 0, P = 0) at large values of the resonance number 
£, because the amplitude of the resonance term due to Eq. (|13D quickly decays when the 
radius, p, decreases. Indeed, at £ = 4 and p = 0.1 the amplitude of the resonance term with 
m = £ = 4 in Eq. is 80 times less than the amplitude of the non-resonant term with 
m = i — 1 = 3. In order to describe the motion in the region near the CGS, we consider the 
initial Hamiltonian (|l]) under the condition X ^ 1. The exact classical equation of motion 
reads. 



Vok^ . 



M 



sin(X - nt). 



(14) 
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FIG. 2. The trajectories in the central resonance cell in the phase space for 6 = and (a) £ = 3, 
e = 5 X 10-"^, (b) ^ = 4, e = 0.05, (c) ^ = 5, e = 0.05. 



Up to the first order in X Eq. (|T^) is. 



—X + u;\l-ecos{nt))X = -^sin(fit). 



(15) 



where e = Vok'^/Mu'^ is the dimensionless perturbation amplitude. If we introduce a new 
dimensionless time, 2r = Qt, then from Eq. (|T3|) we obtain the Mathieu equation with the 
additional right-hand side term in the form. 
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X + a^(l — ecos(2r))X = a£esm(2r), 



(16) 



where = {2/tf. From the theory of Mathieu functionsi it is known that for small e, 
Eq. ([T^) has unstable general solutions at = 1 and = 4 which correspond to the 
resonance numbers, C. = 2 and £ = 1. The additional term in the right-hand side of Eq. 
(p!6|) does not influence the stability of trajectories (see Ref.i, §6.22). At < 1 and small 
enough values of e, the Mathieu equation has periodic solutions which correspond to the 
stable dynamics for resonance numbers £ > 2. In Fig. 2 stable trajectories in the system 
described by the Hamiltonian (|l|) are shown for ^ = 3, 4, 5. Stable region in the vicinity of 
the CGS can be considered as additional "central cells" to those resonance cells shown in 
Figs. 1 (c) - 1 (e). The motion in the central cell is characterized by the following features: 
i) the size of the cell increases with increasing the resonance number, ii) trajectories in 
the cell have the axial symmetry of the order £; iii) as follows from numerical calculations, 
the period of the motion in the central cell along all trajectories in Figs. 2 (a) - 2 (c) is 
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FIG. 3. The trajectories obtained by solution of Eq. (p^ when only the terms up to (a) 
the first order (Eq. (|IB|)) and (b) the fourth order in X are taken into account; £ = 5, 5 = 0, 
e = 0.05. 
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approximately the same (in each figure), and very large. For example, the period of motion 
along the trajectories in Fig. 2 (b) for £ = 4 is (1 — 1.2) x lO^T, where T = 2tt/Q, while the 
period of motion along the trajectories in Fig. 1 (d) is of the order ~ 5 x lO^T. The property 
i) can be described by the fact that the resonance term, which leads to unstable solution 
near the point (X = 0, P = 0), for larger values of i has less infiuence on the dynamics 
because its amplitude at p -C 1 quickly decreases with £ increasing, due to Eq. (|l3l). 

In order to treat the properties ii), iii), we have considered the infiuence of the terms of 
higher order in X in Eq. (|14D on the dynamics described by Mathieu equation (|l^). The 



results for the case i = 5 are shown in Figs. 3 (a), 3 (b). From Figs. 3 (a), 3 (b) it is 
seen that the terms of high order in X change the shape of trajectories, and make them 
symmetrical with the axial symmetry of the order i. This fact follows also from comparison 
of different trajectories in Fig. 2 (c). Indeed, internal trajectories with small values of p 
have the round form, unlike external trajectories with larger values of p, which have the 
form of the pentagon. From comparison of Figs. 2 (c) and 3 (b) it is clear that the terms 
of high order in X (the terms of the order X* where i > i — 1) in Eq.([T^) do not influence 
signiflcantly the dynamics, because the form of trajectories in Fig. 3 (b), computed for the 
approximate model, is practically the same as the form of trajectories in the exact model 
shown in Fig. 2 (c). 

It was established numerically that slowness of the motion in the vicinity of the CGS 
is also result of influence of the terms of high order in X. Thus, the period of motion 
along the external trajectory described by the approximate equation (|16]) in Fig. 3 (a) is 
~ 4 X lO^T. Including into consideration the second order term in X increases the period of 
motion up to ~ 3.3 x lO^T. If we include, as well, the third order term, the period becomes 
~ 3.6 X lO^T. Finally, including the resonance term leads to increase of the period up to 
the value ~ 3.8 x lO^T shown in Figs. 2 (c) and 3 (b). 

Next, we shall analyze the djTiamics in the central cell depending on e. Modiflcation of 
trajectories at increasing the wave amplitude, e, for the resonance number £ = 4 is shown 
in Figs. 4 (a) - 4 (d). Two features in the structure of the trajectories of the central cell 
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can be observed, i) Increase in e shifts the central cell up. ii) The size of the cell increases 
considerably in Figs. 4 (a) - 4 (c) with increasing e in comparison with the case shown in 
Fig. 2 (b), when the value of e is very small. The cell shrinks when e increases, which is 
shown in Fig. 4 (d) for e = 10. Further increase of e destroys the central cell entirely. 
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-0.6 -0.3 




X 



FIG. 4. Trajectories in the central resonance cell at £ = 4, 6 = under influence of the pertur- 
bation with the amplitude, e: a) e = 0.5, b) e = 2, c) e = 5, d) e = 10. 
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Similar features were observed in dependence of the dynamics in the central cell on e for 
the case i — 5, shown in Figs. 5 (a) - 5 (d). Comparison of the data for £ = 5 in Figs. 5 (a) 
- 5 (d) with those for £ = 4 in Figs. 4 (a) - 4 (d) allows us to conclude that the area of the 
central cell increases with increasing i, and chaotization of the motion in the central cell at 
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FIG. 5. Trajectories in the central resonance cell at ^ = 5, (5 = under influence of the pertur- 
bation with the amplitude, e: a) e = 5, b) e = 10, c) e = 15, d) e = 19. 
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larger values of i requires larger values of e. In other words, the motion in the central cell 
becomes more stable with increasing the resonance number, £. In order to understand the 
observed dynamics, we again included into consideration only the first order terms in X in 
Eq. (0), and computed the dynamics for large values of e. The trajectories described by 
the approximate equation for e = 5, £ = 4 and for e = 9.5, £ = 5 are shown in Figs. 6 
(a), 6 (b). The following features can be observed, i) As follows from our calculations. 



2.5 



1.5 
P 

1 

0.5 




(a) 



-1 -0.5 



X 



3.5 



2.5 

P 

1.5 



0.5 



(b) 










1 ^ 








«\;« 
'■I 



0.5 



-1.5 



-0.5 



0.5 



1.5 



X 



FIG. 6. The phase trajectories described by the approximate equation ([T6| ) for 5 = 
and (a) £ = 4, e = 5, (b) £ = 5, e = 9.5. 

the phase portrait shifts up from the point {X = 0, P = 0) under the infiuence of the term 
in the right-hand side of the approximate equation (|16]). ii) Comparison of Fig. 6 (a) with 
Fig. 4 (c) and Fig. 6 (b) with Fig. 5 (b) allows us to conclude that the terms of the high 
order in X in Eq. (0) change the shape of trajectories and restrict the region of stable 
motion, iii) The motion described by the approximate equation (|1^) becomes unstable at 
values of e > where en lies in the interval 2.3 < £3 < 2.4 for £ = 3; 5.5 < £4 < 5.6 for 
£ = 4; and 9.6 < 65 < 9.7 for £ = 5. One can see from Fig. 4 (d) and Figs. 5 (c), 5 (d) that 



the motion described by exact equation (|1J) in the region e > remains stable. Thus, the 
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high order terms in X in Eq. (|T4D stabihze the dynamics at large values of the perturbation 
amplitude, e. 
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(c) ^ = 5, e = 1; (d) i = 5,€ = 5. 
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Let us compare the classical dynamics in the central cell with the dynamics in the other 
cells when the perturbation parameter e is not small. The results of calculation of the 
dynamics in several cells are shown in Figs. 7 (a) - 7 (d). From comparison Fig. 7 (b) with 
Figs. 4 (c), 4 (d) and Fig. 7 (d) with Figs. 5 (a) - 5 (d) one can see that the trajectories in 
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the central cell remain stable, while other nearest cells are completely destroyed by chaos. 
The extremely high stability of trajectories in the central cell can be explained by relatively 
small influence on the dynamics of the terms of high order in X, oscillating with different 
frequencies, because their amplitudes are small at small X (or p). 

III. THE CLASSICAL DYNAMICS NEAR THE CGS IN THE NEAR 

RESONANCE CASE 

Now, let us consider the CGS in the near resonance case, when 5 0. The resonant 
Hamiltonian takes form 

He = i{Suj) + VoJfXp) cos e. (17) 

The stationary points for the dynamics generated by the Hamiltonian ([17D are defined by 
the conditions, 

e = dHe/di = 0, / = -dHe/de = 0. 
Positions of the elliptic stationary points are given by the expressions, 

Vo \J- = T5uJ, = 0, TT, (18) 

ol 

where the sing "— " corresponds to the stable point, with the angle 6^ = with the P-axis, 
and the sing "+" corresponds to the stable point with the angle 6e = it. In the dimensionless 
form Eq. ( pTsD is, 

1 dJM 6 ^^g^ 



pe dp ie 

where pe = kr{Ie). For the positions of the hyperbolic stationary points one has. 



kr{h)] =0, eh = ±|. (20) 



As one can see from Eq. (0), the number of the elliptic stable points in the near resonance 
case, when 5 7^ 0, is finite because the right-hand side of Eq. (|19|) is constant while the left- 
hand side oscillates, and decreases on average. As a consequence, there is a finite number 
of the resonance cells. 
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The motion near the CGS can be described by the approximate equation (|T6D with the 
parameter equal to: = [2/(£ — 5)]^. It is known! that at small e and 6^0 the Mathieu 
equations have the stable solutions at any i including the cases i = 1 and i = 2. 

Let us consider the cases i = 1 and £ = 2 at 5 7^ in detail. In the case i = 1, one may 
use the results of the resonance theory at arbitrary small X and P, because the term of the 
lowest order in X (proportional to X) in the Hamiltonian (|l|) is resonant. Let us suppose 
that the dimensionless radius pe in Eq. (\[% is small, pe ^ 1. Then J[{pe) ~ 1/2, and Eq. 
(H) yields, 

Pe = Te/{26i). (21) 

Thus, shift of the stable elliptic point from the CGS is small, pe ^ 1, when the condition 
e <^ 2\6\i is satisfied. At small values of the wave amplitude, e, shift of the elliptic stable 
point from the point X = 0, P = is proportional to e. One can see from Eq. (|2T|) that at 
i = 1 one elliptic stable point exists at arbitrary small value of e (which follows also from 
the theory of Mathieu functions). 
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FIG. 8. The phase space for the near resonance case, 5 = 0.1, i = I, and (a) e = 0.4, (b) e = 0.7, 
(c) e = 1.2 

When e is small, the phase trajectories are the circles with the center located near the 
CGS. Figuratively speaking, in the case 5 7^ and e ^ 1 there is only one resonant (central) 
cell with an infinite area, because at small enough value of e the equation (|19]) has no other 
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solutions, except for Eq. ([21|), and in the phase space there are no other cells, except for the 
central one. When e increases (we suppose e > 0) and 6 > 0, the stable point shifts down, 
as shown in Fig. 8 (a), because the left-hand side of Eq. (p!9| ) is positive and we should take 
the sign "+" in the right-hand side, which corresponds to the shift in the direction 6 = n. 

When e increases (see Fig. 8 (b), 8 (c)), the dynamical chaos appears, and the area of the 
central cell decreases. As before, we have considered the influence of the high order terms 
in X in the exact equation of motion ([l^) on the dynamics described by the approximate 
equation (plGl) . The approximate equation (|l^) yields unstable solutions when e > ei, where 
|ei| = ^24|5|/5 if 5 > 0, and ei = y^24|5| ii 5 < O.i The parameter 5 = 0.1 yields ei = 0.69.[| 
As one can see from Figs. 8 (b), 8 (c) the central cell remains undestroyed. Thus, the 
nonlinear terms stabilize the dynamics in the near resonance case, similar to the case of 
exact resonance. At e = 1.2 in Fig. 8 (c) one more cell is generated, because the condition 
(pj|) is satisfied for two values of kr. 
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FIG. 9. The phase space for the near resonance case, 6 = 0.1, £ = 2, and (a) e = 0.17, (b) 
e = 0.22 



We also checked these criterion numerically. 
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Unlike the case i = 1, when i = 2 and e is small enough (see Figs. 9 (a)), the stable point 
does not shift from the point {X = 0, P = 0). Instead, in Fig. 9 (b) we observe bifurcation 
at the value e = €2, where €2 can be estimated from the solution of the approximate equation 
(p!6|). Namely, up to the second order in 5, the dynamics becomes unstable at €2 = 25 — 5^/2.1 
Our computed value of €2 lies in the interval 0.185 < £2 < 0.186 which is slightly less than 
the estimated quantity due to the influence of nonlinear in X terms, which are neglected in 
the approximate equation ([16|) . 




FIG. 10. The same as in Figs. 9 (a), 9 (b) but for e = 0.8 

As shown in Fig. 10, at further increase of e, two stable points, formed after bifurcation, 
diverge at larger distance from each other, the chaotic area increases, and additional cells 
appear because the condition (|19|) is satisfied for more number of points (A;r(/e)). 
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IV. STABILITY OF THE QUANTUM GROUND STATE 



Now we consider a stability of the ground state of the quantum harmonic oscillator 
(QGS) under the same conditions as in the classical model. The quantum Hamiltonian is, 

H = ^ + ^-x' + cos(A;x - fit) =Ho + V{x, t), (22) 

where p = —ihd/dx, and the same notation as in Eq. (|l|) where used. Since the Hamiltonian 
(^) is periodic in time, we can use the Floquet theorem, and write the solution of the 
Schrodinger equation in the form, 

ipg^Xjt) = exp{—iEqt/h)uq{x,t), (23) 

where Eg is the quasienergy, ipq{x, t) is the quasienergy (QE) eigenfunction, and the function 
Ug{x,t) is periodic in time, Ug{x,t) = Ug{x,t + T), where T = 2tt/Q. We expand Uq{x,t) in 
the basis of the unperturbed harmonic oscillator, 

oo 

Uq{x,t) = Y.C?,it)Mx), (24) 

n=0 

where coefficients, C^{t), are periodic in time, C^{t) = C^{t + T). Due to periodicity of 
C^(t), the approach based on Floquet states is very convenient for investigation of local- 
ization properties of the quantum system. Namely, if some initial state coincides with the 
quasienergy function localized in some region of the Hilbert space, C„(0) = C^(0), then it 
will remain localized in this region for any time. 

We used the following numerical procedure to calculate the QE states.0Hll The QE 
states are the eigenstates of the evolution operator for one period of the wave field, f/(T). 
In order to build the matrix Unm of the operator U{T) we choose the representation of the 
Hamiltonian Hq. Let us act with the evolution operator on the wave function ip{x, 0), 

U{T)ij{x,0) = ij{x,T), (25) 

and choose the initial state in the form C„(0) = 6n,no- Iii this way we obtain a column in 
the evolution operator matrix, 
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FIG. 11. The plots of averages versus dispersions for h = 0.5, 5 = and (a) £ = 1, e = 0.05, 
(b) ^ = 2, e = 0.05, (c) ^ = 3, e = 5 X 10"^, (d) £ = 4, e = 0.05, (e) ^ = 5, e = 0.05. 
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n,no 



(26) 



where the coefficients, C^''\T), can be obtained by numerical solution of the Schrodinger 
equation (for more detailed discussion see Ref.0). After diagonalization of Un,m, we obtain 
the QE functions, = C^{mT), m = 0, 1, 2, . . . and the quasienergies. Eg. The values 
of matrix elements, Un,m, depend on three dimensionless parameters: the wave amplitude, 
e, the quantum parameter, h = k'^h/Mu, which can be treated as a dimensionless Planck 
constant, and from the ratio Q/u = £ — 6. 
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FIG. 12. The probability distribution for the QGS QE states with the smallest average, Ug, for 
the cases (a) i = 4, (h) £ = 5; e = 0.05, h = 0.5, 6 = 0. 

When 6 = and the amplitude of the wave is small, e ^ 1, most of the QE states are 
divided into almost independent groups, each located in one resonance cell of the Hilbert 
space.i In order to show this let us characterize each QE state by its average, ng = J2n ^l^nl"^, 
and a dispersion, cTg = [J2n{^ ~ )2|C^|Y^^ and plot Hg versus ag (see also Ref.li^l). Such 
plots for different values of the resonance number, i, and for small value of the wave ampli- 
tude, e, are shown in Figs. 11 (a) - 11 (e). The boundaries of the cells are marked by arrows. 
The radius of the external boundary of classical cells in Figs. 1 (a) - 1 (e) corresponds to 
the position of the boundary of the first quantum cell, respectively, in Figs. 11 (a) - 11 (e) 
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with the quantized radius, p„ = \Jlnh. One can see from Figs. 11 (a) - 11 (e) that the QE 
states are mostly located within the quantum resonance cells, because their averages, n^, are 
situated inside the cells, and their widths, cXg, do not exceed the width of the corresponding 
resonance cell. Such states form rows in Figs. 11 (a) - 11 (e). Each cell in the Hilbert space 
in the quasiclassical limit corresponds to 1i classical cells.0 There also exist the QE states 
which do not belong to a particular resonance cell, but rather they belong to the stochastic 
web. These QE states have a "separatrix" structure, i.e. they delocalized over several reso- 
nance cells and have large dispersion, a. Such states are represented by scattered points on 
the diagrams Uq = riq^aq) in Figs. 11 (a) - 11 (e). The structure of these QE states was in 
details discussed in Refs.i'iiii. 




FIG. 13. The same as in Figs. 12 (a), 12 (b) but in the logarithmic scale. 

In this paper, we focus on the QE states which belong to the central resonance cell (QGS 
QE states), because these states are mainly responsible for stability properties of the QGS. 
Such QGS QE states are characterized by small average, n^, and are located in the low part 
of the plots, ng(crg), in Figs. 11 (d), 11 (e). In the cases £ = 1,2 in Figs. 11 (a) - 11 (b) 
these states are absent, which corresponds to unstable dynamics near the CGS in the phase 
classical space, shown, respectively, in Figs. 1 (a) - 1 (b). In the case £ = 3, the area of the 
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stable island in Fig. 2 (a) is much less than the value of the dimensionless Planck constant 
h {h = 0.5). So, in this case, the QGS QE state is absent, too. 

The plots of the probability distribution for the QGS QE states q' with the smallest 
average, nq>, marked in Figs. 11 (d), 11 (e) by arrows, are shown in Figs. 12 (a), 12 (b), for 
the cases i = 4 and i = 5. In the logarithmic scale, these states are illustrated in Figs. 13 
(a), 13 (b). The Husimi functions of these states are shown in Figs. 14 (a) {i = 4) and 14 (b) 
= 5). As one can see from Figs. 12 (a), 12 (b) the QGS QE state is mainly localized in 
the CGS of the harmonic oscillator. One can see from Figs. 12 (a), 12 (b) and 13 (a), 13 (b) 
that the small part of the probability distribution is located at the levels with the numbers 
n = £m, where m = 1, 2, . . .. This can be explained by influence of the resonance terms in 
the quantum equations of motion.! In the resonance approximation the QE states can be 
defined from the set of algebraic equations which in the dimensionless form can be written 
as, 

(fe " ^'') = l(^r.,n+eCl_,, + V;,„_,C^,). (27) 
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FIG. 14. (a) The Husimi functions of the QGS QE state shown in Fig. 12 (a), (b) The Husimi 
functions of the QGS QE state shown in Fig. 12 (b). The cross-sections are plotted from the level 
0.047 with the increment 0.042, h = 0.5, e = 0.05. 
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The matrix elements for n S> 1 can be approximated by the Bessel functions Jm,Q 



n,n+2m+l 



e i 



2 J{n + 1) . . . (n + 2m + 1) 



-J'. 



2m+l 



(28a) 



1 



--J2m{V2nh). 



(28b) 



2 ^{n + 1) . . . (n + 2m) 
As one can see from Eq. (pTf), in the resonance approximation the QE functions have 



the form = with m = 1, 2, . . .. In this case, the particle is allowed to move only 
between the states with the numbers n = im. As shown in Ref.0, a particular form of the 
QE function, = C^^, makes the Husimi functions, illustrated in Figs. 14 (a) and 14 (b), 
symmetric with the axial symmetry of the order i. One can see that the width of the Husimi 
distribution in Figs. 14 (a) and 14 (b), is, AP ~ AX y/h 0.7. 
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Fig. 15 (b) 
h = 0.5. 
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(a) The Husimi functions of the shifted from the ground state QE function, shown in 
The cross-sections are plotted from the level 0.047 with the increment 0.042, e = 5.0, 



^More precise form of matrix elements see in Ref.t 
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Now, let us consider influence of dynamical chaos on the QE states, when e increases. 
There are several QE states localized near the QGS of the harmonic oscillator. Some of 
them are shifted from the level with the number n = 0. They can be associated with shifted 
up classical central resonance cell, when e increases. In this case, the Husimi function in 
Fig. 15 (a) of the QE state with the probability distribution illustrated in Fig. 15 (b) (for 
i — 5 and e = 5) has a similar form to the form of the trajectories in the classical phase 
space shown in Fig. 5 (a). At large enough values of the wave amphtude, (e = 5 in Fig. 15 
(a)), there are no QE states localized in the nearest quantum cells, except for the QE states 
localized in the central cell. This corresponds to the chaotic classical dynamics shown in 
Figs. 7 (b), 7 (d) with the stable island in the center of the phase space. 




O 2 4 e 8 1 o 

wave amplitude 



FIG. 16. The probability to find the system in the QGS of the harmonic oscillator, defined by 
the values jC^^ol^ Q^ state q' mostly localized at the state with n = 0, versus the wave 

amphtude, e, for three values of the effective Planck constant, h; £ = 5. 

The QGS QE states, mostly localized at the CGS of the harmonic oscillator, are of the 
most interest for us, because they mainly define the dynamics of the quantum state initially 
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located at the level with n = 0. The time-evolution of the system with the initial state 
Cn(0) = 6n,no IS defined by the equation, 

CnimT) = Y: Cl*Cl exp{-tE,mT/n), (29) 

where m = 0, 1, 2,.... The amplitude of probability to find the system in the initial state, 
no, is, 

Cn,{mT) = J2 \Cl\^exp{-tEgmT/h). (30) 
g 

Suppose that some QE state with the number q' is mostly localized at the level n = hq, i.e. 
l^ni)P ^ l^noP ^'^^ 1 7^ l'- Then, the term with q = q' dominates in the sum in the 
right-hand side of Eq. (^), and we can write, 

C„„(mT) ^ \C^J^exp{-iE^,mT/h). (31) 

The probability, PnQ^mT), to find the system at the moments tm = mT in the state with 
n = riQ is given by, 

Pr,,{mT) = \G,,{mT)\' ^ \Ci\\ (32) 

The value of P„o (mT) in this approximation is independent of the number of periods passed, 
Pn^irnT) = PnQ. In the next approximation, the neglected terms in Eq. (|30D cause the 
probability P„q slightly oscillate with time. 

In Fig. 16, we present a plot of the probability, Pq = |Cno=ol^) ^"^^ ^^e system in the 
QGS as a function of the wave amplitude, e, if the initial state is the QGS of the harmonic 
oscillator. One can see that the dynamical chaos (the range of large enough e) decreases 
this probability. However, the process of delocalization of the QGS QE state is extremely 
slow when e increases, in comparison with that in other nearest cells. For example, at e = 5 
all QE states in the nearest cells are chaotic (delocalized) which corresponds to the chaotic 
classical dynamics in Fig. 7 (d), while the QE state located at the QGS remains localized 
with the probability Pq ~ 0.56 when h = 0.5, and Pq ~ 0.6 when h = 1.0. From comparison 
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FIG. 17. (a) The QE states mostly localized at the QGS of the harmonic oscillator (QGS QE 
states), h = 0.5, (a) e = 0.05, (b) e = 0.5, (c) e = 5.5 (d) e = 6.0, (e) e = 7.6, (f) e = 9.2. The 
boundaries of the quantum cell are marked by arrows. 
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of different curves, for h=0.1, h = 0.5 and h = 1.0, one can note the following features, i) At 
small values of e, increase of h leads, on average, to decrease of stability of the QGS. ii) The 
QGS, at large values of /i {h = 1), is more stable under the influence of chaos (the range of 
large enough e) than that at small values of h {h = 0.1). We should note, that oscillations 
of Po{mT) in time should increase when Po{mT) decreases. This happens in the region of 
large enough e in Fig. 16, because in this case the influence of neglected terms in Eq. ( |30D 
becomes more signiflcant. 

More information about the stability of the QGS can be extracted from the analysis of 
the structure of QE states located at the QGS at different values of the wave amplitude, e, 
shown in Fig. 17 (a) - 17 (f) (for h = 0.5). The QGS QE state shown in Fig. 17 (a) for 
small value of e (e = 0.05) is similar to the separatrix QE states,! because it has the regular 
structure (compare, for example, with Fig. 17 (c)), and its maxima are located near the 
quantum separatrices, indicated in Fig. 17 (a) by arrows. Thus, the QGS QE state also 
po=^s=e= the properties of the ".eparatrix' QE states, considered in Retl. 

The separatrix QE states are of the quantum natureB because they are delocalized over 
several resonance cells. These QE states provide the tunneling between the cells when 
chaotic regions in the phase space are negligible small, and the classical particle can not 
practically penetrate from one resonance cell to another. The separatrix QE function mostly 
localized in the QGS is a particular one, and it is different from other separatrix states 
studied before. On the one hand, it is delocalized over several resonance cells (see Fig. 17 
(a)) as other separatrix QE functions. On the other hand, this particular QE function is 
mostly concentrated on the QGS of the harmonic oscillator, unlike the other separatrix QE 
functions. These "contradictory" features of the QGS QE state deflne the dynamics: on the 
one hand, the system mainly remains localized in the QGS, but on the other hand, a small 
part of the probability distribution can tunnel to the states with large n, located in the other 
resonant cells. When we increase the parameter h, the separatrix QE states become more 
delocalized, and the probability to tunnel to other cells increases, which explain decrease of 
Pq with increasing h in Fig. 16 (compare the different curves in the region of small e). 
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At intermediate values of e, when 1 < e < 7, a stability of the QGS increases with 
increasing h. However, we can not increase h indefinitely because when h becomes larger 
than the resonance area in the phase space, the QGS becomes unstable. Thus, at = 5, 
£ = 4, e = 2 (see the classical phase space in Fig. 4 (b)) and at/i = 5,£ = 5,e = 5 (Fig. 5 
(a)) no localized QGS was found. 

An increase of Pq with increasing the wave amplitude, e, when e is small (e < 0.5 for 
h = 0.5 and e < 1.0 for h = 1.0), shown in Fig. 16, is a consequence of a partial localization 
of the separatrix QE state (see Fig. 17 (b)) under the influence of chaos explored in Ref.0 
In this case, the QGS QE function, shown in Fig. 17 (b), loses its "separatrix" features. 
Further increase of e, makes the QGS QE state more delocalized. However, as one can see 
from Figs. 17 (c) (e = 5.5, h = 0.5) and 17 (d) (e = 6.0, h = 0.5) delocalization takes place 
mainly over the nearest oscillator states with small numbers, n. At e > 5.5 {h = 0.5) the 
oscillations appear in the dependence Pq = -Po(e), in Fig- 16. Thus, the QE function at 
e = 6 in Fig. 17 (c) is less localized than the QE function at e = 7.6 shown in Fig. 17 (d). 
At large values of e (e > 7.6), practically all QE states are delocalized, which is the quantum 
manifestation of chaotization of the classical central cell in the phase space. 
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FIG. 18. Time-evolution of the dispersion a = cr{m), where m = t/T, for three values of the 

wave amplitudes: (a) e = 0.05, (b) e = 0.5, (c) e = 7.6, and for h = 0.5, 1 = 5,5 = 0. 



In order to illustrate a non-monotonic character of localization of the QGS as a function 
of the wave amplitude at small e, we computed the dynamics of the quantum state initially 
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concentrated on the ground state of the harmonic oscillator, C„(0) = Sn,o, using Eq 
Time-evolution of the dispersion, 



a{mT) = J2 |C„(mT)|2(n - n(mT))2, (33) 
y n 

where n{mT) = J2n \Cn{TTT'T)\'^n is the average, is presented in Figs. 18 (a) - 18 (c), for three 
values of e. When the wave amplitude, e, is small (Fig. 18 (a)), a small part of the wave 
packet can propagate to large values of n due to effect of diffusion via the separatrices as 
shown in Fig. 19 (a). Similar tunneling effect of the wave packet between the resonance cells 
via the separatrices was explored in Ref .1 In spite of a small probability of tunneling to other 
cells, contribution of this part to the dispersion, a, is significant because it is proportional 
to {n — n)^, where n — n ^ 1. At e = 0.5 the separatrix QE states are destroyed by chaos, 
as shown in Fig. 17 (b), and QGS becomes more localized (see Fig. 19 (b)). This leads to a 
significant decrease of a in Fig. 18 (b) in comparison with the case of small e, shown in Fig. 
18 (a). Further increase of e, up to the value e = 7.6, results in delocalization of the QGS, 
as shown in Fig. 19 (c), and the dispersion, cr = a{mT), in Fig. 18 (c) becomes large again. 
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FIG. 19. Probability distribution in the system with initial condition (7^(0) = dn,no and (a) 
e = 0.05, rriT = 3 X lO^T; here |CoP = 0.746, (b) e = 0.5, mT = lO^T; |CoP = 0.982, (c) 
e = 7.6, mT = 5000T; |CoP = 0.137, |Ci|^ = 0.423. The boundaries of quantum cells (quantum 
separatrices) are marked by arrows; h = 0.5, i = 5, 5 = 0. 

Comparison of Figs. 19 (a) and 19 (c) allows us to conclude that delocalization of the 
QGS at very small, and at large values of e is of a different nature. In the former case, the 
diffusion is caused by the separatrix QE states. These states are the quantum objects because 
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they are delocalized, and provide tunneling between the resonant cells-Q As a consequence 
of the quantum nature of the separatrix QE states, increase of the dimensionless Planck 
constant, /i, leads to delocalization of the separatrix QE states and decrease of localization 
of QGS at small e, shown in Fig. 16. On the other hand, when e is large we observe 
delocalization caused by chaos, which is manifested in the irregular form of the probability 
distribution, shown in Fig. 19 (c). 

As follows from Fig. 16, in order to make the QGS more stable at small e one should 
decrease the Planck constant, h. A plot of the dispersion, a, versus time for h = 0.1 is 
presented in Fig. 20 (a). The system remains in the ground state with the probability 
Pq = 0.996. However the dispersion is still large enough because the particle can propagate 
with small probability to the levels with large n ^ 1, due to the diffusion via the separatrices 
which can be seen from the plot of the probability distribution presented in Fig. 20 (b). 

Another way to increase the stability of the ground state at small e, is to destroy the 
separatrix QE functions by choosing the non-resonant value of the wave frequency, so that 
6 = i — Q/uj 7^ 0. For a non-resonant case {6 = 0.01), the plot of the dispersion as a function 
of time is presented in Fig. 21 (a), and the probability distribution at time m = t/T = 3 x 10^ 
is illustrated in Fig. 21 (b). One can see from Fig 21 (a) that introducing a detuning, 5 7^ 0, 
results in considerable improvement of the stability of the QGS in comparison with the case 
of the exact resonance, illustrated in Fig. 18 (a). Thus, in order to make the ground state 
more stable at small values of e one must detune the system from the exact resonance. 

The minimal value of detuning, 6, required to destroy the separatrix structure, and to 
make the ground state more stable, can be estimated from the quantum equations of motions 
in the resonance approximation (]27|). The term, proportional to 6n destroys the separatrix 
QE states as shown in Ref.i This term becomes significant when it becomes of the order 
of the relation, e/h. On the other hand, the separatrix QE function must, at least, occupy 
the two separatrices. Thus, we can estimate the number, n = ns, of the oscillator's state at 
which the separatrix structure is destroyed, namely ns = e/6h. The separatrix QE functions 
decay exponentially with increasing n in the region n > n^.i For the parameters in Fig. 21, 
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ns — 50, which is much less than the position of the first separatrix rii — 385 shown in Fig. 
20 (b) by arrow. So, the separatrix QE states at these parameters do not exist, and we do 
not observe tunnehng from the QGS to other resonance cells. 
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FIG. 20. (a) Time-evolution of the dispersion, a = cr(m), where m = t/T, and (b) probability 
distribution at time m = 3 x 10^ for small value of h, 6 = 0, h = 0.1, e = 0.05 
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FIG. 21. (a) Time-evolution of the dispersion a = a{m), where m = t/T, and (b) probability 
distribution at time m = 3 x 10^ for the near resonance case, 6 = 0.01; h = 0.1, e = 0.05 
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Decreasing the value of e, in the near resonance case, makes the QGS more stable. Unlike 
the near resonance case, in the case of the exact resonance, a stability of the QGS at e <^ 1 
is independent of the wave amplitude, because in this case the localization properties of 
the QGS are defined by the structure of separatrix QE function, which in the resonance 
approximation (when e is small) is independent of e.i One can see from a comparison of Fig. 
21 (a) with Fig. 20 (a) {h = 0.1, e = 0.05) that the dispersion in the near resonance case 
(Fig. 21 (a)) is much less than that in the exact resonance case (Fig. 20 (a)), in spite of the 
probability to stay in the ground state, Pq, for these two cases does not differ significantly 
(Po = 0.99897 at 6 = 0.01 and Pq = 0.996 at 5 = 0). The reason is that the dynamics in the 
exact resonance case is mainly determined by the separatrix QE function, and in this case 
there is a small probability for particle to tunnel to the high oscillator levels with n ^ 1, as 
shown in Fig. 20 (b). 
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FIG. 22. The same as in Fig. 16, but for the near resonance case, 5 = 0.01. 

When the wave amplitude e increases, the condition 6^0 becomes less significant. 
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This can be explained by less influence of the term proportional to S on the dynamics in 
comparison with influence of the wave with the amplitude e in the region where the value 
of kr is relatively small (see the classical Hamiltonian in Eqs. (^), and quantum 

Eq. (|2^)). In Fig. 22 we plot the function Pq = Po(e) for the near resonance case. In Figs. 
23 (a) - 23 (c) we compare the results for the exact resonance case {6 = 0) with those for 
the near resonance case, when 6 = 0.01 and 6 = 0.1. One can see from these flgures that 
the dynamics in the vicinity of the QGS in the near resonance case is similar to that in the 
exact resonance case, except for the region of small e, which was discussed above. 




wave amplitude wave amplitude wave amplitude 

FIG. 23. Plot of Pq versus e for the exact resonance case, when 6 = (solid line and open 
circles), and two curves for the near resonance cases: 5 = 0.01 (dashed line and filled squares), 
5 = 0.1 (dot-dashed line and crosses); i = 5, (a) h = 1, (b) h = 0.5, (c) h = 0.1. 

In the quantum case, similar to the classical one, at very small e and flnite S, there always 
exists the QGS QE state at any i including the cases i = 1 and £ = 2. In Fig. 24 we present 
a plot Pq = Po(f) the cases i = 1 and i = 2 when 6 = 0.1 and h = 0.5. As one can 
see from Fig. 24 the QGS QE state in these two cases becomes unstable at considerably 
less values of e than that in the case £ = 5 in Fig. 23, which corresponds to the classical 
dynamics in Figs. 8-10. In the case i = 1 the stable point shifts down from the region 
X = 0, P = in Figs. 8 (a)- 8 (c), and in the case i = 2 the stable point at e ~ 0.2 becomes 
unstable (see Figs. 9 (a), 9 (b)) and 10). The quantum manifestation of this process is a 
rapid decay of the value of Pq in the region e > 0.2 in Fig. 24. 
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In conclusion, the classical dynamics in the vicinity of the point {x — 0, p — 0) in 
the classical phase space and quantum dynamics in the vicinity of the ground state of the 
harmonic oscillator in a field of a monochromatic wave, is explored. Both resonance and 
near resonance cases are analyzed. It is shown that at small e and finite detuning from the 
resonance, 5, the quantum ground state is always stable. In the case 5 — 0, the dynamics 
is unstable for the resonance numbers i = 1,2. Stability of the classical dynamics in the 
central cell and stability of the quantum dynamics near the ground state of the harmonic 
oscillator under the influence of chaos is analyzed. It is shown, that under certain conditions 
{£ > 2, 5 = 0, e ^ 1, h ^ 1) the presence of chaos makes the quantum ground state more 
localized. Increase of the quantum parameter, h, at intermediate values of e, enhances 
localization of the QGS considerably. Experimental observation of discussed results in the 
system of an ion trapped in a linear ion trap and interacting with laser fields represents a 
significant interest for understanding the stability properties of this system. 




O 0.2 0.4 0.6 0.8 1 

wave amplitude 

FIG. 24. Plot Pq = -Fb(e) for the resonance numbers £ = 1 (open circles and solid line) and 1 = 2 

(filled squares and dashed line) in the near resonance case, 6 = 0.1; h = 0.5. 
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